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We determine both the in-plane and out-of-plane dynamics of viscoelastic membranes separating 
two viscous fluids in order to understand microrheological studies of such membranes. We demon- 
strate the general viscoelastic signatures in the dynamics of shear, bending, and compression modes. 
We also find a screening of the otherwise two-dimensional character of the response to point forces 
due to the presence of solvent. Finally, we show that there is a linear, hydrodynamic coupling 
between the in-plane compression modes of the membrane and the out-of-plane bending modes in 
the case where the membrane separates two different fluids or environments. 



I. INTRODUCTION 

Membranes and biopolymers form many of the most 
basic structures of plant and animal cells. These fun- 
damental building blocks frequently occur together in 
complex structures. In animal cells, for instance, the 
outer membrane is often strongly associated with a net- 
work of filamentous actin, one of the most prevalent 
proteins in the cell. This actin cortex is viscoelastic 
and it contributes significantly to the response of whole 
cells to external stress. Many prior physical studies, 
both theoretical as well as experimental, have concerned 
the structure and dynamics of simple membranes Q and 
interfaces Much less is known about complexes of 
membranes with biopolymers. Recent experiments have 
demonstrated the ability to both construct and probe in 
vitro models consisting of lipid membranes with attached 
actin biopolymerQ. By a microrhcology technique, the 
material properties of these micrometer-scale viscoelastic 
films could be measured. Here, we calculate the dynam- 
ics, and thereby the response, of viscoelastic membranes. 
We demonstrate, among other things, the general signa- 
tures of viscoelasticity in the dynamics of both shear and 
bend. We also find a coupling of bend and compression 
modes when there are different environments on the two 
sides of the membrane, as is the case for the surface of a 
biological cell. These effects have important implications 
for previous and ongoing microrheological studies of both 
the model biopolymer-membrane complexes, as well as 
real cells. 

Microrhcology 1^, ||, ^ studies the rheological proper- 
ties of a material by the use of small probe particles, 
which can either be actively manipulated by external 
forces[0, or imaged while subject to thermal fluctua- 
tions. It is possible, for instance, to extract the vis- 
coelastic moduli from observations of the fluctuating po- 
sition of a small (Brownian) particle embedded in the 
medium [|[ ^, ||, |[ 0. The technique holds out great 
promise as a new biological probe measuring the mate- 
rial properties of living cellsH. In effect, the judicious 



application of this technique may permit the creation 
of a "rheological microscope," providing new insights 
into the regulation and time-evolution of the mechani- 
cal properties of various intracellular structures over the 
course of the cellular lifecycle, or in response to vari- 
ous external stimuli. Such studies of the cyctoplasm are 
already underway ||ll|, [l^ . One essential feature of the 
cell is the cell membrane, an essentially two-dimensional 
lipid bilayer incorporating a wide variety of dissolved 
proteins and anchored to a cytoskeletal network. Mi- 
crorheological studies of artificial biopolymer-membrane 
complexes Q have been performed, and efforts are under 
way to extend such techniques to the membranes of real 
cells In order to study the rheology of the cell mem- 
brane using microrheological techniques, the previously- 
studied methods used to extract rheological measure- 
ments from thermal fluctuations in three-dimensional 
samples [|[ |[ ^ |l^, |l^ need to be extended to prob- 
lem of a viscoelastic membrane coupled to a viscous sol- 
vent. 

In this paper we consider such an extension of these 
ideas, with an eye toward not only cellular microrhcol- 
ogy, but also the investigation of a wide variety of sys- 
tems in soft physics wherein a viscoelastic membrane is 
coupled to a viscous fluid (typically water). Examples 
include emulsions, vesicles, and Langmuir monolayers. 
We calculate the position response to a force of a small 
rigid particle embedded in the complex, soft, viscoelastic 
medium. Using the Fluctuation-Dissipation theorem p6| 
we can then compute the autocorrelations of the probe 
particle's position - a quantity accessible in experiments. 
As in previous work in this field 0, |, ||, 0, ||, careful 
attention must be paid to the full linear mode structure 
of the system. The observed thermal fluctuations of the 
probe particle are in response to all such modes that cou- 
ple to the particle position thus the full mode spectrum 
must be determined in order to interpret microrheological 
data. This is in distinction to more traditional rheology, 
which is a linear response measurement of the system 
to applied shear strain and therefore measures the shear 
modulus directly. In particular, in this membrane study 
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we are required to calculate the hydrodynamic flows (in 
the low Reynolds number limit) in the viscous fluid above 
(superphase) and below (subphase) the membrane gen- 
erated by membrane distortions. The understanding of 
these flows is essential to calculating the modes of the 
combined system and, hence, the response function of 
the probe particle. This point has been previously rec- 
ognized by Brochard et al. ^ in their calculation of the 
decay of membrane bending perturbations when coupled 
to a viscous fluid, and was also used in the interpretation 
of the recent experiments on actin-coated membranes . 

Qualitatively, we find that the introduction of an 
viscous subphase and/or superphase introduces a new 
length scale over which shear waves in the two- 
dimensional membrane decay due to the viscous damp- 
ing of the three-dimensional fluid. The appearance of 
a new length scale is, of course, not surprising given 
that the ratio of a two-dimension membrane shear mod- 
ulus and a three dimensional fluid shear modulus 
G{uo) — —iusri yields a length. In fact this length has 
been commented upon previously in the context of vis- 
cous films in contact with a viscous solvent In 
the case of membrane dynamics and microrheology, it 
sets a probe particle-dependent, high frequency hmit, be- 
yond which the dynamics are controlled entirely by the 
solvent^. In general, this length sets a limiting range of 
the two-dimensional strain or velocity field. 

In addition to the role of membrane-liquid coupling in 
the shear modes, we perform a similar analysis of bend- 
ing modes of the membrane. We reproduce the Brochard 
et al. results for the mode structure in the case that 
subphase and superphase fluids have the same viscosity, 
and we assess the out-of-plane response function in the 
manner outlined briefly above. More remarkably, we find 
that linear independence of the in-plane membrane com- 
pression modes and the out-of-plane bending modes is 
destroyed when the z — > —z symmetry of the problem 
is broken by a membrane (lying in the xy plane) that 
separates two fluids of different viscosities. We go on 
to calculate the response function in this more general 
case as well as in the symmetric problem. The asymmet- 
ric problem has clear physical relevance to dynamics of 
cellular membranes as well as to vesicles and microemul- 
sions. Perhaps the most dramatic example of the broken 
symmetry occurs in the study of Langmuir monolayer dy- 
namics where the subphase (typically aqueous) and the 
superphase (air) have viscosities which differ by many 
orders of magnitude. 

Finally we point out that it is well-known that in mem- 
branes whose equilibrium shapes are not flat, there is a 
linear coupling of bending modes to in-plane distortions 
due only to the nontrivial geometry of the surface. We 
do not discuss this sort of coupling in the present work 
and restrict the present analysis to the dynamics of flat 
membranes leaving the role of curvature in microrheology 
to a later paper. 

The remainder of this paper is organized as follows: 
In section II we discuss the coupled modes of the mem- 



brane viscous fluid system. We use this analysis of the 



mode structure in section HI to determine response func- 
tion of a rigid particle embedded in the membrane in 
the symmetric case where in-plane and out-of-plane mo- 
tions decouple. We then break the subphase/superphase 
symmetry in section IV and reanalyze the response func- 
tion for the case where out-of-plane forces can generate 
in-plane motion. Finally, in section ^ we compute the 
predicted position autocorrelations of a particle embed- 
ded in a membrane whose material properties are typical 
of the various classes of systems mentioned above. We 
pay particular attention to the system of an actin-coated 
lipid bilayer that has been investigated by the Strasbourg 
group 1^. We conclude and propose further experimental 
and theoretical work in VI. 



II. THE MODES OF THE SYSTEM 

We now determine the modes of the membrane cou- 
pled hydrodynamically to the (typically aqueous) sub- 
phase and superphase. The essential calculation of these 
involves solving the equations of motion for the fluid 
above and below the membrane given a certain defor- 
mation mode of that membrane. We shall do this for 
the three independent modes of membrane deformation. 
These modes include two in-plane modes, shear and com- 
pression, as well as one out-of-plane, bending mode of 
the membrane. We consider first the modes of in-plane 
membrane deformation, i.e. those that do not involve 
curvature of the membrane. As long as the membrane is 
fiat in equilibrium, these modes are all linearly indepen- 
dent. We shall see, however, that in systems having fiuids 
of different viscosities in the subphase and superphase, 
there is a linear-order, hydrodynamic coupling between 
in-plane compression mode and the out-of-plane mem- 
brane bending. However, if the two fluids are identical 
restoring the z — > —z symmetry of the problem, the 
three modes mentioned above are all linearly indepen- 
dent. We will study the modes of the broken-symmetry 
problem in a later section. 

The membrane lies in the xy plane as shown in flgure 
|l|. The strain field u is a two-dimensional vector lying in 
the plane of membrane. We now need to determine the 
fiuid velocity field above the membrane {z > 0) associ- 
ated with this shear wave. Working in the limit of zero 
Reynolds number we solve the Stokes equation 



VP, 



(1) 



where the three dimensional vector v is the fluid velocity 
field and P is the hydrostatic pressure which enforces the 
fluid incomprcssibility 



V-v = 0. 



(2) 



These equations must be solved subject to the bound- 
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ary conditions 

v(x,z = 0,t) = ^u(x,i) (3) 
lim v(x, 2,i) = (4) 

z >oo 

reflecting the stick boundary conditions of the fluid at 
the surface of the membrane and requirement that the 
fluid velocity field go to zero at large distances from the 
membrane. 

A. Shear deformation 

We choose membrane coordinates so that the shear 
wave propagates in the x direction and the deformation 
is in the y direction. Thus the simple, in-plane shear 
deformation of membrane is described by the strain field 

u(j/,t)=yC/oe^(«--"*) (5) 

From the symmetry of the problem we look for a 
solution of Eq. ^) with boundary conditions given by 
Eqs. (|,||)of the form 

V = -zc.C/oy/(z)e*(«^-'^*\ (6) 

where j{z) is an unknown function satisfying the condi- 
tions: /(O) = 1, lim^ ,oof{z) = so that the ansatz, 

Eq. (j^) satisfies the requisite boundary conditions. The 
Stokes equation demands that the vorticity of fluid flow, 
V X V, satisfies Laplace's equation 

(V X v) = 0. (7) 




FIG. 1: The flat membrane considered in this paper. As 
shown in (I) the membrane lies in the xy plane of out coordi- 
nate system. The in-plane displacement field u is defined in 
the plane of the membrane. The fluid superphase z > and 
subphase z < are not shown. In (II) the edge-on view of the 
membrane shows the vertical displacement h of the membrane 
(solid line) from its fiat, equilibrium shape (dotted line). 



Using our ansatz, we find that unknown function f{z) 
satisfies the differential equation 

§-,'1 = 0. (8) 

Along with the boundary conditions given above for f{z), 
we find a solution for the velocity field in the region above 
the membrane {z > 0) 

v(x, z, t) = -icjC/oye*«^-l«l^e-*'^*. (9) 

The fluid velocity in the subphase (z < 0) is similar 
with z — > —z. Returning to the Stokes equation we 
find that there is no pressure gradient associated with 
this fluid motion. We see that the shear flow induced in 
the three dimensional viscous liquid phases decays over 
a distance comparable to the wavelength of the in-plane 
shear mode. This viscous damping by the surrounding 
fluid introduces a new decay length for shear waves in 
the membrane even in the case that the membrane were 
perfectly elastic, i.e. n{oj) = /xq, a real constant. More 
on this later. We now turn the other in-plane mode of 
the membrane, the longitudinal or compression mode. 

B. Compression mode 

We now apply a longitudinal compression wave in the 
membrane having a strain held 

u(x, t) = xUo cos{qx)e-"^* . (10) 

and determine the associated fluid velocity field in the 
superphase (z > 0). 

From an examination of the fluid flow near the 
fluid/membrane boundary, we note that the compression 
mode injects a sinusoidally varying vorticity field that 
is directed along the y-axis and varying in the x direc- 
tion. Since the vorticity must satisfy the Laplace equa- 
tion [Eq. (^] and since the boundary conditions require 
it to vary sinusoidally in the x direction, the fluid vortic- 
ity must decay exponentially into the fluid, i.e. in the z 
direction. Based on these considerations, we expect that 
the vorticity O = V x v takes the form 

f7 = T-iycos(ga;)e"l«l^ (11) 

in the superphase. The fluid velocity must vanish at large 
distances from the membrane so we have chosen a decay- 
ing exponential in the z direction in Eq. (|l^). The con- 
stant T, with dimensions of time is, as yet, undetermined. 
It will be selected to enforce the stick boundary condi- 
tions of the fluid at the membrane's surface. From the 
above equation and incompressibility we flnd the differ- 
ential equation obeyed by the z-component of the fluid 
velocity field in the superphase 

{dl + d^,)v, = ^e-\'^\'sm{qx). (12) 
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Once again, we use the ansatz for v^: ^^(x, ^)e*"* = 
f{z) sin(qx) and find a differential equation of for f{z) 
of the form 



(13) 



The boundary conditions at the membrane and at infinity 

require that /(O) = Um^ >oo f{z) = 0. 

The homogeneous solution of the above boundary con- 
dition vanishes upon the application of the boundary con- 
ditions on f{z). This leaves only the particular solution 
to give the solution for Vz in the superphase. Integrating 
the fluid inconipressibility equation, Eq. (^, then gives 
the accompanying solution for Vx- We find the fluid ve- 
locity to be given by: 

Vzix,z,t) = -icjC/oe"'"* [1 - |g|z]e-l«l^cos(ga;)(14) 



Vz{x,z,t) = -iujUoqe "^*ze '*'^sin(gx). 



(15) 



Returning to the Stokes equation we can calculate the 
pressure field associated with the above fiow. It is inter- 
esting to note that there is a sinusoidally varying pres- 
sure field at the surface of the membrane. We find that, 
for the compressional wave in the membrane introduced 
above, the pressure at the upper surface of the membrane 
(z — > 0+) takes the form 



P{x,t) = 2ir]LL!qUosm{qx)e 



(16) 
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It appears from the hydrodynamic flow field shown in fig- 
ure |, from and pressure field calculated above in Eq. (|l^ 
that there should be a membrane deformation along its 
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FIG. 2: The membrane compression mode with associated 
flow field in the superphase. The membrane, seen edge-on 
as a dotted hne undergoes a compression wave of wavelength 
27r. The vector field represents the fluid flow in the superphase 
and shows the exponentially decaying vorticity. 



normal accompanying the membrane compression mode. 
Of course, in the symmetric case where the fluids in the 
superphase and subphase are identical, such a normal 
deformation must vanish by symmetry. As soon as this 
symmetry is broken, however, there should be a mem- 
brane height fluctuation in response to the longitudinal 
modes of the membrane. This is, in fact, incorrect. A 
complete calculation of the zz component of the fluid 
stress tensor at the surface of the membrane shows that 
the pressure term is exactly canceled by the viscous stress 
arising from the gradient of the upward fluid velocity. 
Nevertheless, such a linear order, purely hydrodynamic 
coupling between membrane compression and bending 
modes exists and is mediated by the surrounding fluid as 
we shall see in our analysis of the bending modes. There 
the flows generated by the membrane bending generates 
a viscous shear stress at the surface of the membrane 
which couples to the in-plane compression mode. 

This effect is unrelated to the coupling of shear and 
bending modes in membranes with finite mean curva- 
ture [? ]. The new linear coupling could be realized in 
a variety of experimental situation. Among these are 
microemulsions, polymersomes, cell bodies, and Lang- 
muir monolayers. These systems can break the symme- 
try of the problem even in the limit of flat membranes 
(or membranes with a large radius of curvature leading 
to a vanishingly small geometric coupling between shear 
and membrane bending) by having fluids of differing vis- 
cosities on either of the membrane. We return to this 
point later; we point out only the fundamental asymme- 
try of the coupling. Bending modes generate in-plane 
compression, but in-plane compression does not produce 
membrane bending. 



C. Bending mode 

In this section we recapitulate the results of Brochard 
et al. on the effect of hydrodynamics on the dynamics 
of membrane bending modes. To do this we apply a 
sinusoidal height fluctuation to the membrane of the form 



h{x,t) = ;i,e*(«^-'^*) 



(17) 



where h{x,t) measures the displacement of the membrane 
surface in the direction normal to its surface, i.e. the z- 
direction. See flgure [l[ Once again, using the Stokes 
equation we calculate the fluid flow in the superphase 
generated by such a membrane displacement. The vortic- 
ity generated by the membrane motion must still satisfy 
Laplace's equation the symmetry of the problem admits 
flows only in the xz plane. We assume that the solution 
takes the form 



v(x, z, t) = [vxiz)x + vz{z)z] e'(9-"*). 



(18) 



Using the Laplace equation for the vorticity and the 
incompressibility condition we flnd two differential equa- 



tions for the unknown functions Vx{z),Vz{z) 

dvz _ 
dz 



n dVrr 



dz^ 



iq- 



dz-^ 



dz 







(19) 
(20) 



Combining Eqs. ( p9| , po| ) in order to eliminate Vx we arrive 
at a differential equation governing Vz alone 



d^Vz 



,<Pvz_ 

dz^ 



q^Vz = 0, 



(21) 



dz4 

which has the solution: 

Vziz) = Ciel«l^ + C2e-I«l^ + Cg^el"!^ + dze'^'^^'. (22) 

To satisfy the boundary condition at infinity we set Ci = 
C3 = 0. After integrating Eq. ( [l^ ) to obtain a solution for 
Vx{z) and applying the stick boundary conditions at the 
surface of the membrane, we return to Eq. ( [T^ ) to write 
the solution for the fluid flow field in the superphase 



Vzix,Z,t) = 



(23) 



-iujho [1 + \q\z] e-l9ke^(9^-'^*' (24) 

(25) 

These flows are sketched in figure |3[ As required, the 
component of the fluid velocity field tangent to the mem- 
brane vanishes. 

The pressure gradient associated with the above flow 
field can be calculated directly from the Stokes equation. 
We find, upon setting the pressure to zero at infinity, that 
the pressure field at the surface of the membrane takes 
the form 



P{x, z = 0, i) = 2itj|g|7y/ioe^(«^-"'*'. 



(26) 



We note in passing that fact that the amplitude of the 
pressure oscillation is proportional to q combined with 
the restoring force on the bent membrane depending on 
wavevector as q* leads to the well-known result that the 
decay rate of bending modes increases as the third power 
of wavevector. 



III. THE RESPONSE FUNCTION 

We now turn the calculation of the response function. 
The calculation is similar in spirit to calculation of the re- 
sponse function in three dimensions. We wish to model 
the response to a force of a particle embedded in the 
membrane at the origin. Given the level of description of 
the membrane, i.e. a continuum, single-component, two- 
dimensional viscoelastic medium, there is no distinction 
between the application of a force to the membrane itself 
over an area the size of the particle, and the applica- 
tion of the same force to a rigid particle embedded in the 
medium. Because of our reliance on a simple, coarse- 
grained description of the membrane, we necessarily ne- 
glect any new hydrodynamic modes of the membrane 



that may, in fact, be present in the physical system. For 
example, if we were to consider a two-component mem- 
brane composed of lipid and an elastic protein network 
anchored in the lipid bilayer, there should be, in anal- 
ogy to the three-dimensional gel system, a "free drain- 
ing" mode associated with the diffusive relaxation of net- 
work density. Nevertheless, we expect that such modes 
will be irrelevant at higher frequencies where our single- 
component description of a membrane should be valid for 
many systems of experimental interest. Details related to 
the more complex descriptions of the membrane internal 
structure and their effect on microrheological measure- 
ments we defer to a future publication. In addition we 
do not consider the effect local perturbations of the mem- 
brane structure (and consequently its viscoelastic proper- 
ties) due to the introduction of the probe particles. It has 
been shown experimentally ||l^, |l3| and theoretically |l4| 
in three dimensional systems that such local perturba- 
tions can be important in the one particle response func- 
tion, but that interparticle response functions do not de- 
pend on such effects to leading order. 

We apply a force on the membrane that is localized 
at the origin of the coordinate system and is sinusoidally 
varying in time. We then compute the response of each 
Fourier mode of the membrane deformation along with 
the associated fluid motion. By integrating over these 
motions we determine the motion of the point at the 
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FIG. 3; The membrane bending mode with associated flow 
field in the superphase. The membrane, seen edge-on as a 
solid line undergoes a bending wave of wavelength 2ti. The 
undeformed membrane is shown edge-on as a dotted line. The 
vector field represents the fluid flow in the superphase. The 
component of the fluid velocity field tangent to the mem- 
brane generates a shear stress that couples to the in-plane 
compression mode. This shear stress changes sign across the 
membrane, however, only the flow in the superphase is shown 
above. 
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origin of the membrane. The response of each Fourier 
mode in the viscoelastic membrane is determined from 
force balance 

= [fJid'^Ua + (m + \)dadf)Ui3\ - Si^nd^h + crjjz=o + ft- 

(27) 



In the above equation, the Greek indices run over the 
coordinates in the plane of the membrane, x, y, while the 
Latin indices run over all three coordinates. The first 
term in square brackets on the RHS of Eq. (27) repre- 



sents the in-plane membrane viscoelastic force per area 
due to shear and compression of the membrane. The 
quantity d^^ projects out the membrane coordinates. The 
membrane viscoelasticity is described by two frequency- 
dependent Lame constants consistent with a isotropic 
continuum. These terms should in fact be expressed as 
integrals of over the strain history of membrane. We will, 
however, suppress the frequency dependence of the Lame 
constants (and thus the viscoelasticity of the membrane) 
until we rewrite the force-balance in the frequency do- 
main. The force per area associated with the out-of-plane 
displacement of the membrane, h(xa,t), is given by the 
second term on the RHS of the above equation, which 
represents the restoring force on the membrane due to its 
bending rigidity, k. Once again, in a viscoelastic mem- 
brane, we may assume that k is frequency dependent, 
but we suppress this dependence for now. 

Finally the last term on the RHS of Eq. (|^) is the 
force per area exerted on the membrane by the external 
fluid of the sub- and superphases. The fluid stress tensor 
takes the well-known for an incompressible, Newtonian 
fluid 



(28) 



The fluid velocity and hydrostatic pressure fields that ac- 
company any membrane motion has been determined in 
the previous section so we may rewrite the fluid stress 
solely in terms of the membrane displacement fields: 
and h. Specializing to the components the fluid stress 
tensor required in Eq. (|2^) and noting the stick bound- 
ary conditions at the surface of the membrane we may 
simplify the form of the full fluid stress tensor to 



f 



(29) 



Because we are still considering symmetric case, the a 
force acting in the plane of the membrane excites only 
the in-plane deformation modes discussed above: com- 
pression and shear. Similarly, the component of applied 
force acting along the membrane normal only generates 
bending deformations. Using this decoupling of the in- 
plane and out-plane of plane modes we may calculate the 
response to in-plane forces and forces along the mem- 
brane normal (z) separately. We begin with in-plane 
forces 



A. In plane response 

We rewrite the force balance equation, Eq. ( p7| ) spe- 
cializing to the case of an in-plane force, /q(x) and 
we Fourier transform in two-dimensional membrane sub- 
space. We arrive at 



fiq'^Ua{c[,uj) + {fJ. + X)qaqi3Up{q,u}) - ^Waiq^uj) 



/a(q,w), 
(30) 



where we have defined = dzVa\z=o- Projecting out the 
longitudinal and transverse parts of the above equation 
and writing — UaQa, u 
<ja<ji3, we have 



Pipup with Plf^ = S„p - 



1 T 



f{q)aqa 



(31) 
(32) 



Using the results of section II, in which we have com- 
puted the fluid flows associated with the longitudinal and 
transverse modes of the membrane, we can determine the 
form of w in terms of the u. In this way we can write 
Eqs. (^l|),(|3^) solely in terms of u and thus solve for the 
membrane displacement in terms of the externally ap- 
plied force, f . From the previous section one finds: 



Waqa 



— ia;w'"(q) 
iwu^(q)|q|. 



(33) 
(34) 



Combining the above equations with the force balance 
equations, Eqs. (^,(^) we solve for u(q) = qqu^{(\} + 

{5 —qq)-v7 (q) . We now integrate over the modes of the 
system to determine the displacement of the point at the 
origin in response to the applied force at that point. The 
rigidity of the tracer particle at the origin is accounted 
for by cutting off the wavevector integral at the inverse 
radius of the embedded particle, |<?max| ~ 1/a or equiva- 
lently 



/„(q,c^)-Fo„e-*'^*e((?„,ax-|q|), 



(35) 



where Fq is a constant vector, which for definiteness we 
take to be in the x direction. We choose an order one 
numerical prefactor in this relation between Qmax and 1/a 
so that the response function reproduces the standard 
stokes drag on a spherical particle in the limit that the 
membrane elasticity vanishes. Having done this we find 
that the position response of the embedded tracer sphere 
(which is along the x by symmetry) takes the form 



U:c(x,w) = Fq 



<fq 



(2^)2 ^(2;. + A)g2_Ifl|q| 
1 — cos^ (h 



(36) 



n - viqi J 

where the angle (j) is defined hy f ■ x = cos (p. Clearly 
the integral is the response function that we seek and 



7 



via the Fluctuation-Dissipation theorem, it contains the 
information necessary to determined the experimentaUy 
measured power spectrum. Performing the integrals we 
arrive at the final form of the in-plane response function. 
Since it is diagonal in the in-plane indices we suppress 
them and write: 



a||(w) = 



1 
+ 



log 1 



(37) 



2// + A 



log 1 + i 



.2(2/i + A) 



We having chosen the wavevector cutoff here to be gmax — 
2/(37ra) so that in the limit of no membrane elasticity: 
fi, A — > with /i/A finite, the response function reduces 
to the standard Stokes drag result from low-Reynolds 
number hydrodynamics. 



B. Out— plane— response 

We now analyze the out-of-plane motion in a similar 
way by first introducing a Fourier component of the bend- 
ing deformation of the membrane of the for m of E q. (17). 



Using our hydrodynamic results of section II C 
bination with Eq. ( ^7|) and computing the hydrodynami- 
cally induced stress using Eq. (ES) we determine the am- 



plitude of a bending mode of the membrane in response 
to an arbitrary force. We then specialize to the case of 
a point force and, integrating over the available bending 
modes of the system, we determine the out-of-plane re- 
sponse of the point at the origin to a force on it directed 
along the membrane normal. 

We make use of the fluid flow field associated with the 
bending deformation of the membrane as discussed in 
From this solution and Eq. 



section 



[IC 



3|) we find that 

at the surface of the membrane there is a non-vanishing 
component of the fluid stress tensor which takes the form 



a,, — 2iujh, 



i{qx—Lot) 



(38) 



It should be pointed out that although there is a nonva- 
nishing fluid stress component, cr^^, this stress changes 
sign across the membrane and exactly cancels in the sym- 
metric case. We return to this point later in the paper 
when we examine the asymmetric problem in detail. Fur- 
thermore, we point out that if we were to consider a mem- 
brane as having finite thickness and thus being able to 
support internal strains in which the deformation gradi- 
ent normal to the membrane's surface, i.e. a membrane 
that can support a stress strain of the form u^z, then the 
bending mode would hydrodynamically couple at linear 
order to such internal deformations of the membrane. We 
do not pursue this point in the current paper. For our 
present purposes, it is enough to note that a force normal 
to the plane of the undeformed membrane couples only 
to the bending modes in the symmetric case. 

From the z component of the force balance equation 
[Eq. (|27|)] we find that a force in the z direction with a 



sinusoidal dependence on Xa and t generates a sinusoidal 
membrane bending deformation with amplitude 



h{q) 



Kq'^ — 2iujri\q\ 



(39) 



Integrating over such forces in q-space, determine the re- 
sponse to a point force at the origin using the wavevector 
cutoff introduced in Eq. (|35|)we find the out-of-plane re- 
sponse function 



2k 



dp- 



1 



i6 



(40) 



where the dimensionless parameter S = 27Tr'^a^uj'r]/{4:K) 
is the ratio of the viscous stress urj/a to the membrane 
bending stress k/c^ at the length scale of the probe par- 
ticle. We have now completed our discussion of the re- 
sponse function of a bead embedded in the membrane for 
the symmetric case. In this case the complete motion of 
the particle, x(ci;), in the membrane is a linear combina- 
tion of the in-plane motion due to in-plane components 
of the applied force and out-of-plane motion due to the 
component of the externally applied force normal to the 
membrane. Thus the most general solution of the me- 
chanical problem to linear order in membrane displace- 
ments takes the form 

Xt{uj) = a\i{uj){Sij - ZiZj)fj{uj) + a{uj)z!z{t^) (41) 

where the response functions a_L{ijj) and az{uj) are given 
by Eqs. (^^ and ( ^0| ) respectively and is the exter- 
nally applied force responsible for the motion. We will 
later use the Fluctuation-Dissipation theorem to explic- 
itly compute the implications of this result for the ex- 
perimentally observed position fluctuations of the probe 
particle. First, we turn to the case of the asymmetric 
problem to examine the coupling of in-plane membrane 
compression modes to the out-of-plane membrane bend- 
ing modes and the implications of such coupling on the 
response function. 



IV. ASYMMETRIC SYSTEM 

Here we examine the system in which the membrane 
separates fluids of differing viscosities. In order to study 
membrane dynamics in this asymmetric case we return to 
the question of the hydrodynamic stresses on the mem- 
brane with a particular emphasis on the coupling of com- 
pression and bending (out of plane motion) modes me- 
diated by the hydrodynamic coupling through the sub- 
phase/superphase. From our previous calculation of the 
fluid velocity fields associated with the various deforma- 
tion modes of the membrane, we determine the viscous 
shear stress on the membrane by demanding the conti- 
nuity of the shear stress Uaz, a — x,y sX the membrane- 
fluid boundary. In the fluid we note that this component 
of the stress tensor takes the form 



r]dzVa 



(42) 



8 



We note that daV^ = 0,a = x,y due to the fluid bound- 
ary conditions on the membrane. Due to the hnearity of 
the hydrodynamics, we can combine the solutions of the 
flow flelds from the longitudinal (compression) and bend- 
ing modes to write the full shear stress associated with 
a linear combination of those membrane deformations. 
We may deal with the transverse membrane modes sep- 
arately since they do not couple to the bending degrees 
of freedom as does the compression mode. Finally, we 
note that the shear flows on the opposite sides of the 
membrane generated by the longitudinal in-plane modes 
contribute additively, while the flows coming from the 
out-of-plane deformation of the membrane produce shear 
stresses on either side of the membrane that tend to can- 
cel. Both of these shear stresses, however, are propor- 
tional to the viscosity of the fluids on either side of the 
membrane; thus there will be a residual shear stress act- 
ing on the membrane when there are fluids of differing 
viscosities on either side of it. In this case, the shear 
stress associated with the compression mode will be pro- 
portional to the sum of the viscosities on either side of 
the membrane, while the shear stress associated with the 
bending mode will be proportional to the difference in 
those viscosities. Taking these observations into account, 
we find that the shear stress from the fluid at the surface 
of the membrane may be written as 



we arrive at 



(43) 



In a similar way we may compute the normal stresses 
on the membrane due to the hydrodynamic flows. The 
stress component normal to the membrane is 



(44) 



where P is the hydrostatic pressure computed from the 
fluid velocity field and the Stokes equation. Here we find 
that only the bending mode generates a nonvanishing 
normal stress component. The compression mode pro- 
duces a pressure variation across the surface of the mem- 
brane that exactly cancels the z-derivative of the vertical 
velocity field. We write the normal stress as 



lz=0 



2iL0 (S??) \q\h{q). 



(45) 



Returning to the Fourier-transformed force balance 
equation we write out the components in the q z sub- 
space. The z equation is identical to Eq. (39). We use 



this result to eliminate the dependence in the q equation, 
which takes the form: 



Bq^u^{q)- 2quj 



1 



isign((?) (S?7) u'^iq) - -{A7^)h{q) 



(46) 



where f (q) is the externally applied force and B = 2fi + X. 
Solving for the longitudinal part of the displacement field 



u'^iq) 



Bq^ - 2iuj\q\ {Y.ri) 

c^|g|(Ar/)/,(g) 

[Kg4 - 2iijj (Sry)] [Bq^ - 2iuj\q\ (E??)] ' 



(47) 



Combining the above result for the longitudinal part of 
the in-plane displacement field with the previously cal- 
culated transverse part of the in-plane displacement field 
as well as the out-of-plane perpendicular displacement, 
we can write the trajectory of any point on the membrane 



'R.(Xa,ijj) — zh{Xa,Oj) + Vi{Xa,Uj), (48) 

which in the Fourier transformed variable takes the form 
Ri{qa,uj) = Zih{qa,uj) + qiU^{qa,uj) + {S^p - qiqp) up{qa,uj). 



(49) 



We now calculate the one- and two-particle response 
functions for the asymmetric membrane using the above 
results. To do this, we localize the applied force at the 
origin by taking f{xa, ^ij) = exp(— jwi)f(5(x) and we cal- 
culate the displacement field of a particle attached the 
membrane at some other location X using the relation 



R(;f,c.)= j -f±^n{q^,u)e^^-^ . 



(50) 



From this calculation we determine the response function 
tensor aij{X,uj) defined by the equation 



(51) 



so that aij {X , oj) measures the displacement of a point at 
X on the undeformed membrane in response to a force 
applied at the origin of the coordinate system. As be- 
fore, we take into account the finite size of the particle 
at the origin (where the force is applied) and the tracer 
particle at X by cutting off the wavevector integrals at a 
maximum wavevector proportional to the inverse particle 
radius. We choose the coefficient of proportionality to be 
the same as we have already used above. For concrete- 
ness we take the displacement vector of the observation 
point at X to be along the i-axis. It is, of course, triv- 
ial to rewrite the resulting expressions in an arbitrary 
reference frame. 

Using Eqs. (|47|),(^), our previous solutions for the 
bending and transverse, in-plane membrane deforma- 
tions in Eq. (|5^) we first determine the displacement of 
the tracer particle in x direction. The resulting angular 
integrals are simply written as Bessel functions and we 
find 



R^{X,uj) 



,iX,u)U+a,,iX,iu)f, (52) 
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where the components of the response tensor are given 
by the integrals over the magnitude of the wavevcctor q 



Oixz{X,Uj) 



1 

4^ JO 



''^ \ Uq\X\) - JMX\) ,,,,, 



+ 



t Bq^ - 2iuj\q\{i:ri) 
J2{q\X\r 



liq^ - ioj\q\{i:ri) 

iuj{Ar]) 1''^^°' 



277 



dqq 



Uq\x\) 



{Bq^ - 2iw|g|(Er])] 
(54) 



[k^^ _ 2ia;|g|(Sr?)]' 

The most interesting point of this result is the presence of 
a finite Uxz in this problem. We sec that, due to the hy- 
drodynamic coupling through the subphase/superphase, 
a force applied normal to the membrane at the origin 
produces in plane displacements at a different point X . 
This in-plane motion is due to the coupling of membrane 
compression modes to the bending of the membrane gen- 
erated by the applied force in the ^-direction. The ap- 
pearance of this coupling is clearly proportional to the 
viscosity difference between the subphase and the super- 
phase and vanishes in the symmetric membrane prob- 
lem. This hydrodynamic coupling is only observed in the 
two point measurement; in the limiting case where one 
observes the displacement of the particle at the origin due 
to a force applied to it (the \X\ — > limit of the two- 
point measurement which is clearly identical to the one 
point measurement in this calculation) the off-diagonal 
part of the response function vanishes since Ji(0) = 0. 
The vanishing of the off-diagonal component is clearly 
required by the azimuthal symmetry about the z-axis of 
the single particle problem. 

The remaining nonvanishing components of the re- 
sponse function tensor are easily calculated in an analo- 
gous manner. We find that these components are given 
by 



{X,w) = — 



Mq\x\) 



dqq 



+ 



Mq\X\) + Uq\X\) 
Bq^ - 2iw|g|(Ery) 

Uq\x\) 



^5) 



^iq^ - iw\q\{Y.r]) 



for in-plane motion due to an in-plane force perpendic- 
ular to the separation vector of the two particles. For 
displacements normal to the plane of the membrane we 
find 



a. 



.(A',w) 



z/ a 



dqq 



M<i\x\) 



nq'^ - 2i(j{Y,r])\q\ 



(56) 



It should be noted that for large values of their argu- 
ments, Jo (a;) + J2 (a;) is dominated by Jo (.t) — J2 (a;) , so the 
response function along the line of centers between the 
two particles is dominated by the compression modulus 
while the response function for in-plane motion perpen- 
dicular to the line of centers is controlled by the shear 



modulus as one would expect. In the next section we 
turn to the fluctuation dissipation theorem to compute 
the experimentally observable correlated thermal fluctu- 
ations of two tracer particles embedded in the membrane. 



V. POSITION CORRELATION SPECTRA 

Having computed the response function to an applied 
force for both single particle and two-particle systems, 
we now turn to the question of the what is the experi- 
mentally accessible quantity. There are two basic types 
of microrheological measurements that are possible. In 
active microrhcology in which the response function is 
directly probes via the linear response measurement of 
the displacement of a tracer due to a force applied either 
to that particle (one particle measurements) or to an- 
other particle embedded in the membrane (two particle 
measurements). The predicted response functions mea- 
sured in these experiments have been directly calculated 
in this paper. It is, however, also possible to use the 
correlated thermal fluctuations of a single particle or two 
particles embedded in the membrane to access the same 
rheological information via the Fluctuation Dissipation 
theorem. In this section we concentrate on the correlated 
motion of two particles embedded in the membrane in or- 
der to highlight the effect of the hydrodynamic coupling 
between the bending and compression modes that, we 
believe constitute a novel effect. The quantity that we 
compute is 

Sij{X,Lj) = j {Ri{X,t)Rj{0,0))e''^*dt. (57) 

By the Fluctuation-Dissipation theorem this correla- 
tion function can be written in terms of the imaginary 
part of the response functions calculated above as 



Sij {X , LV) 



<.(A',a;), 



(58) 



where 0!"j{X, w) is the imaginary part of the response (in 
the i*^^ direction) of a particle embedded in the membrane 
at A" to a force in the j"^ direction. It only remains for us 
to compute the remaining integrals over the magnitude 
of the wavevector q to determine the correlation spectra. 

For concreteness we put the vector defining the separa- 
tion of the particles along the x direction. We first com- 
pute the correlations of the in-plane motion of the parti- 
cles perpendicular to their line of centers and along their 
line of centers by calculating Syy{X, oj) and Sxx{X, cj) re- 
spectively. In all of the following calculations we assume 
that the in-plane shear modulus is dominated by the real 
part and that this is frequency independent. In other 
words we consider the membrane to act like a perfectly 
elastic sheet. Other calculations can, of course, be per- 
formed with the formulae given above. We present these 
results merely as an example of the correlation functions 
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for this particularly simple case. More complex assump- 
tions about the viscoelasticity of the membrane, presum- 
ably based on microscopic models of the membrane, can 
of course be incorporated in the formalism provided and 
the resulting integrals can be performed. 

We may write the imaginary part of the response func- 
tion in the following form for the motion perpendicular 
to the line of centers as 



1 



47r /iT 



■\X\la 



d z 



2J+(z) 

z2 + t'"^ 



(59) 



where we have defined the functions J^{z) to be the fol- 
lowing combinations of Bessel functions of the first kind: 



J±(Z) = Jo(z)± J2(Z) 



(60) 



and T and r' are a frequency-dependent functions of the 
form 



2uj{J:t])\X\ 
B 



(61) 
(62) 



Using the same definitions, we write the imaginary part 
of the response function for motion along the line of cen- 
ters as 



1 



c\X\/a 



J+{z) f^i\^ 2J-{z) 



B) z'2 + t'-^ 

(63) 



The physical interpretation of r is clear: it measures 
the ratio of the distance separating the two particles to 
the screening length of shear modes in the membrane 
due to the coupling to the viscous subphase and super- 
phase. The other function, r', has an analogous meaning 
in terms of the compression modes of the membrane. Up 
to an overall scale factor the shape of the correlation spec- 
trum obeys a distance-frequency scaling relation so that 
spectra obtained at different interparticle separations can 
be collapsed unto a master curve by rescaling distances so 
they are measured in terms of the frequency-dependent 
screening length. The breakdown of this scaling relation 
can be used as a diagnostic of the appearance some sig- 
nificant frequency-dependence of the membrane elastic 
moduli. 

The only distinction between the two response func- 
tions shown above is interchange of the roles of the func- 
tions J^{z). The effect of this switch is that the compres- 
sion response dominates the long-length scale response 
along the line of centers while the shear response of the 
membrane controls the long-length scale response per- 
pendicular to the line of centers. This observation if in- 
tuitively obvious. We plot the resulting correlation spec- 
tra as a function of the dimensionless variable r in the 
limit that the compression modulus, B = 2fj, + X is much 
larger than the shear modulus in the membrane and that 



the particle size is much smaller than the interparticle 
separation, \X\ ^ a. The plots of the undimensionalized 
correlation functions: lij = Anfi^ Sij / {2kBT{'Eri)\X\) are 
shown in figure ^ 

We now consider the same calculation for the corre- 
lated out-of-plane fluctuations of the particles. Once 
again using the Fluctuation-Dissipation theorem with 
the appropriate component of the response tensor, 



1 



a. 



c\X\/a 



dz- 



^6 _j_ ^11- 



(64) 



where we have defined a new dimensionless, frequency- 
dependent variable in analogy to r and r' above. In 
this case, t" measures whether the viscous damping of 
the bending mode on the length scale of the interparticle 
separation is relevant at the frequency w. It is defined by 



2uj{T.ri)\X\- 



(65) 



We can write the undimensionalized correlation spectrum 
4^ ^i:K^\X\^S,^:,{X,u})/{2kBT{Y.rf}) in terms of an inte- 
gral that is performed numerically. The resulting correla- 
tion spectrum is shown in figure ra where the independent 
variable is r". It may be noted the frequency, separation 
scaling property of the fluctuations of the particles both 
along their line of centers and perpendicular to their line 
of centers fails for the case of the vertical motion. This 
breakdown of the scaling reflects the elementary result 
that bending energy is harmonic in the Laplacian of the 
vertical displacement rather than in the gradient of the 
displacement field as in the case of in-plane deformation. 
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FIG. 4: The graph of the dimensionless interparticle correla- 
tion function (see text) for motion along the line of centers 
Ixx and motion perpendicular to the line of centers lyy as 
a function of the dimensionless variable r in the combined 
limits of small particle size compared in the interparticle sep- 
aration and high membrane bulk modulus compared to the 
membrane shear modulus. 
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Finally we consider the correlations between the out- 
of-plane motion at one point and in-plane motion at an- 
other point. These correlations do not exist in the sym- 
metric membrane which has a vanishing oixz component 
of the response tensor. However, due to the existence of 
such a term in the asymmetric case, these unexpected 
correlations will be present in the thermal fluctuations of 
the compressible membrane. 

Writing the imaginary part of the a^z component of 
the response tensor. 



(Ar?)u;|A'| 



c\X\la 



d zJi{z) 



(66) 



where we have used the definition of t' from Eq. ( |6^ ) 
and have introduced the dimensionless, interparticle- 
separation dependent quantity 



(67) 



B\X\- 



which measures the relative importance of the bending 
modulus and the compression modulus at the length scale 
of the interparticle separation \X\. 

As we have done before, we plot the undimcn- 
sionalized correlation spectrum versus the dimension- 
less frequency variable, r'. In this case the dimen- 
sionless correlation variable takes the form: Ixz = 
-S^z{\X\,u})t:B'^ /{kYiT{/^'q)). We see immediately that 
there are now a family of curves generated by changing 
the dimensionless parameter k' . In figure |^ we plot two 
representative curves from this family: One at a low value 
of K,' (10~^) and one at a high value of this parameter, 
k' = 10^. In addition, due to the interparticle separa- 
tion dependence of this parameter, there is, once again. 



no time-distance scaling as was found for the in-plane 
motion. This correlation function is qualitatively differ- 
ent from the other discussed above in that, at a fixed 
interparticle separation, it changes sign as a function of 
frequency. For concreteness we momentarily assume that 
the superphase is less viscous than the subphase, then, 
at low enough frequencies the vertical motion of one par- 
ticle is correlated with the displacement of the other one 
along their line of centers. This is to be expected from 
the hydrodynamic flow flelds plotted earlier. The upward 
motion of one particle leads to an increase in the separa- 
tion of the two particles. At higher frequencies, however, 
-2t\he size of the rolls in the fluid changes so that the mo- 
lion becomes anti-correlated. Lastly, it is important to 
note that the overall magnitude of the x z correlations 
depends on the in-plane bulk modulus. In highly incom- 
pressible membranes, these correlations vanish as XjB^ . 

We now turn to one final power spectrum calculation. 
Unlike the ones discussed above, we now consider ex- 
plicitly a viscoelastic membrane. Based on the work of 
the Strasbourg group, we take as an example, an actin- 
coated lipid membrane whose viscoelastic properties are 
dominated by the actin coat. We further specialize to the 
high-frequency limit where the actin rheology is domi- 
nated by single chain dynamics. In this frequency range 
the effective shear modulus of the membrane takes the 
form: 



= ^a(-i'i)^/^ 



(68) 



where /ia is a (real) modulus scale and tD = lot is a dimen- 
sionless frequency. We now use Fluctuation-Dissipation 
theorem in conjunction with the in-plane response func- 
tion given in Eq. ( |37| ) in which we set the membrane shear 
modulus to that given above in Eq. (|6^) and the mem- 
brane compression modulus to infinity. This latter choice 
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FIG. 5: The graph of the dimensionless interparticle correla- 
tion function (see text) for motion perpendicular to the mem- 
brane, Izz as a function of the dimensionless variable r" in 
the limit of small particle size compared in the interparticle 
separation. 



FIG. 6: The graph of the dimensionless interparticle correla- 
tion function (see text) for motion perpendicular to the mem- 
brane, /a;z as a function of the dimensionless variable r" in 
the limit of small particle size compared in the interparticle 
separation. 
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is based on the near-incompressibility of the underlying 
hpid bilayer to which the actin network is attached. 

After some algebra we determine that the single par- 
ticle power spectrum is given by 



where the function H{uj) is given by: 



(69) 



Hiu;) 



'37r\ 
cos I — arctan 



cos 



-^sin (7^ ) log 



^l + sin(f)/?M; 
l + /32(cj) + 2/3(w)sin( ^ 



The function H{uj) depends on frequency only through 
the dimensionless quantity 



L 



-1/4 



(71) 



In the limit of high frequency (small /3) it is simple to 
show that the power spectrum scales with frequency as 
Lu~'^ as is expected for simple Brownian motion. One can 
show that in this limit the power spectrum takes the form 



kBT 



as a; ■ 



00. 



(72) 



On the other hand, at intermediate frequencies high 
enough so that Eq. ( |68| ) accurately describes the mem- 
brane rheology, but low enough so that /? is small, H is 
independent of frequency (up to logarithmic corrections) 
and we find that power spectrum the tracer particle po- 
sition fluctuations decays as a different power law with 
frequency. In this frequency range we find that 



-7/4 



(73) 



We also compute the power spectrum of out-of-plane 
fluctuations for the membrane. This result does not de- 
pend on the frequency-dependent shear modulus of the 
actin. Rather, it depends on only the bending modulus 
of the membrane, k, which, in principle, is also a com- 
plex, frequency-dependent quantity. In particular, for 
the actin-coated membrane, the bending modulus takes 
into account the viscoelastic result of the actin in ad- 
dition to the usual bending energy of the lipid bilayer 
so one should expect this quantity to have a complex 
frequency response. Inasmuch as this more microscopic 
mechanical issue has not been satisfactorily resolved, we 
will simply assume that the bending modulus is a real, 
frequency-independent quantity. 

The remaining calculation, of course, is done anal- 
ogously to the calculation presented above. We use 
Eq. ( |40| ) and the Fluctuation-Dissipation theorem to 
compute the power spectrum of z-fluct nations. To do so 
we need to consider the imaginary part of the response 
function by performing the integral in Eq. (^). This in- 
tegral depends on the dimensionless parameter 5 which is 



linearly proportional to frequency. In the low frequency 
limit, we find that to leading order the imaginary part of 
integral takes the form: 



Im[a^(Lj)] = 



2k 



(74) 



From the Fluctuation-Dissipation theorem we then find 
that the power spectrum in the low frequency limit takes 
the form: 



(70) 



1 /e) 2 \ 1/3 



(75) 



Note that the above result is independent of tracer par- 
ticle size. The power-law decay of fluctuations with the 
above exponent has already been calculated by Zilman 
and GraneklHl 

In the high frequency limit, the effect on the tracer of 
the membrane bending stress on the dynamics is domi- 
nated by the viscous stress coming from the fluid. The 
power spectrum for Brownian fluctuations of the bead 
reduces to that of a free Brownian sphere, 



fceT 



as u> 



(76) 



The full power spectrum of the out-of-plane tracer 
particle fluctuations is shown in figure]^; this plot demon- 
strates the cross-over from the low-frequency, bending 
stiffness dominated dynamics, to the high frequency dy- 
namical regime controlled by the viscous stresses in the 
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FIG. 7: The out of plane tracer particle fluctuation spectrum 
in arbitrary units plotted against a dimensionless frequency. 
See the deflnition of S immediately following Eq. (^) . At low 
frequencies, where the dynamics of the tracer is dominated by 
the membrane bending stiH'ness, the —5/3 power law obtains. 
At high frequencies where the fluid viscosity dominates the 
dynamical response of the particle , this power law crosses 
over to the —2 exponent, which is expected for the Brownian 
motion of a free particle. 



13 



VI. SUMMARY 

In this paper we have examined the dynamics of flat, 
viscoclastic membranes either immersed a viscous, New- 
tonian fluid, or separating two Newtonian fluids of differ- 
ing viscosities. We have paid particular attention in this 
analysis to dynamical issues related to microrheological 
measurements. Thus, we have calculated in some detail 
the correlated fluctuation spectrum of two rigid particle 
embedded in the membrane. In addition we have com- 
puted the response function of a single particle in such 
a membrane. Both of these calculations can be directly 
applied to analysis of experimental data, whether it be 
from a actual linear response measurement of the force- 
position response of a tracer particle (e.g. with magnetic 
particles, or laser tweezers) or from the measurement of 
the autocorrelations or interparticle correlations of tracer 
particles undergoing thermal, Brownian motion on the 
membrane. 

The range of application systems is quite large. At 
the moment we only restrict the equilibrium shape of the 
viscoelastic membrane to be flat, although we expect the 
results in this paper to be applicable to curved systems 
as long as the particle motion and interparticle separa- 
tion in the case of the two-particle measurements remains 
small compared to the radius of curvature. The above 
work can be applied to lamellar phases in aqueous sur- 
factant systems, lipid bilayers, emulsion droplets, "poly- 
mcrsomcs"(as long as they are large enough - see the 
above comments regarding curvature), Langmuir mono- 
layers, and cell membranes. 

One of the principle results of this work is the pre- 
diction of a novel cross-correlation between the vertical 
motion of one particle and the in-plane motion of an- 
other one embedded in the same membrane. Such cou- 
pling between the out-of-plane motion of one and the 
in plane motion of the other particle is mediated purely 
by the hydrodynamic interaction through the surround- 
ing liquid and vanishes (by symmetry) if the Newtonian 
fluids above and below the membrane are identical. Such 
cross-correlations should be observable in a large variety 
of systems since in a number of general multiphase and in 
particular biological systems membranes separate fluids 
of different viscosities. This effect, for example, should 



be present in studies of cell membranes. It should be 
most dramatic, however, in cases of high membrane com- 
pressibility and large fluid viscosity differences across the 
membrane. The best system then to look for such cross- 
correlations will be a compressibility Langmuir mono- 
layer. 

A number of extensions of this work are both possi- 
ble and interesting to pursue. The flrst is to study the 
role of equilibrium membrane curvature on the dynamics. 
Since it is already known that there is a purely geomet- 
ric coupling between in-plane strains and out-of-plane 
bending modes, one should find a rich structure in the 
cross-correlations of in-plane and out-of-plane motion 
resulting from the interplay of the geometric and hydro- 
dynamic coupling of these modes. In addition, such an 
analysis is necessary to extend the methods of microrhe- 
ology to highly curved surfaces. Other important ex- 
tensions of the present theory include the analysis of the 
coupling of a viscoelastic membrane to a viscoelastic bulk 
material. The cell membrane coupled to the viscoelastic 
cyctoplasm is an important realization of such a system. 

Along the lines of applying these results to cellular mi- 
crorheology, we have to acknowledge that we have studied 
a one-component, continuum model of the membrane, 
whereas the physical cell membrane is a highly hetero- 
geneous material on the submicron length scale. It is 
therefore important to examine a more complex, multi- 
component model for the membrane. Such more detailed 
dynamical models have been previously studied in the 
context of three-dimensional microrheology. These stud- 
ies in according with very general arguments have shown 
that the there are necessarily extra dynamical modes in 
the heterogeneous models. However, based on the three- 
dimensional work, we expect that the principle effect of 
these extra modes will be to change the correlation spec- 
tra only at the lowest measurable frequencies. 
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